Phase transitions of the g-state Potts model on multiply-laced Sierpinski gaskets 
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We present an exact solution of the g-state Potts model on a class of generalized Sierpinski fractal 
lattices. The model is shown to possess an ordered phase at low temperatures and a continuous 
transition to the high temperature disordered phase at any q > 1. Multicriticality is observed in 
the presence of a symmetry-breaking field. Exact renormalization group analysis yields the phase 
diagram of the model and a complete set of critical exponents at various transitions. 
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I. INTRODUCTION 

The study of models on hierarchically organized lat- 
tices has enriched our understanding of critical phenom- 
ena at phase transitions [U-IH. I n many cases, the renor- 
malization group (RG) transformation of the partition 
function can be carried out exactly [2, H, 0] ■ Exhaustive 
search of the RG fixed points and a complete characteri- 
zation of the flow in the Hamiltonian space then become 
possible. Various conjectures on the universality of the 
critical exponents and their dependence on dimension- 
ality and other model parameters can be tested. The 
intuition gained from the study can be used to guide the 
interpretation of numerical results and experiments. 

In this work, we analyze an exact solution of the q- 
state Potts model [8J on a class of fractal lattices with and 
without a symmetry-breaking field. These lattices, first 
introduced by Menezes and MagalhasesQ, are simple ex- 
tensions of the well-known Sierpinski gasket. They can 
be constructed iteratively following the scheme shown in 
Fig. 1. The Potts model (including the Ising model at 
q = 2) on the original Sierpinski gasket does not order ex- 
cept at zero temperature @, Il0l - tl2| . This may seem sur- 
prising at first sight since the gasket has a fractal dimen- 
sion d( = In 3/ In 2 = 1.58 greater than the lower critical 
dimension d\ c = 1 from field theory. The puzzle was_ 
tially resolved by Gefen, Mandelbrot and Aharony 
who noticed that, in addition to the fractal dimension, 
other topological features of a lattice also affect the ex- 
istence of a finite temperature transition and its critical 
properties. 

The laced Sierpinski gasket can be considered as a hy- 
brid of the Sierpinski gasket and the Berker lattice, which 
does allow a finite temperature transition with a Potts 
Hamiltonian 0, 0]. As in the case of the Berker lattice, 
the coordination number of corner sites on the lattice 
grows with iteration. Presence of such vertices gives rise 
to a hierarchical set of "integration points" that medi- 
ate order among subsystems at low temperatures despite 
their tenuous contact. We construct the phase diagram 
and calculate the critical properties based on exact RG 




FIG. 1. (Color online.) Iterative construction of the multiply- 
laced Sierpinski gasket at m = 2 in two steps and the renor- 
malization group transformation. Potts spins at the corners 
of the triangle (open circles) are used to define the restricted 
partition functions at a given order k. Solid circles indicate 
the internal vertices to be integrated out to yield the renor- 
malized interactions. 



transformations derived from previous work by Borjan et 
al.[12|. Interestingly, when a symmetry-breaking field is 
introduced, we observe multicritical behavior and a line 
of continuous transition between two disordered phases. 

The paper is organized as follows. The fractal lattice 
and the Potts model is introduced in Sec. 2. In Sec. 3 we 
present the iterative relations for the partition function 
and discuss topological features of the RG trajectories 
under the mapping, including the relevant fixed points. 
Section 4 contains calculations of the phase diagram and 
the critical exponents. A brief summary is presented in 
Sec. 5. 



II. POTTS MODEL ON LACED SIERPINSKI 
GASKET 

Figure 1 shows the iterative construction of the fractal 
lattice in two steps. In the first step, the triangle at a 
given order k is triplicated and linked together at the 
three {k + l)th order internal vertices (solid circles). In 
the second step, the structure is replicated m-fold and 
then laced at the three (fc + l)th order corner vertices 
(open circles) to form a (k + l)th order triangle. The 
zeroth order triangle is chosen to be a stack of m unit 
triangles laced at the three corner sites. Each pair of 
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sites are ro-tuply connected. 

Let Nj~, Bk and T k be the total number of sites, bonds 
and basic triangles on a fcth order triangle, respectively. 
From the above definition, we have N k +i = 3mNk — 6m+ 
3, -Bfc+i = 3mB k and T k +i = 3mT k , with Nq = 3, Bo = 
3m and Tq = m. Hence 





(3m) 



fc+i 



3m — 1 



+ 2,B k = (3m) fc+1 ,T fe - m(3m) fe . 



(1) 

Since the linear size of the triangle doubles upon each 
iteration, the fractal dimension of the lattice is given by 
di = ln(3m)/ln2. 

The coordination number or the degree of a fcth order 
corner vertex is given by C k = 2m fe+1 . The coordina- 
tion number of an internal vertex at the same order is 
D k = 2Ck-i = Am k . On a fractal lattice of order n, 
there are three nth order corner vertices. The number 
of fcth order internal vertices can be calculated from the 
iterative relation N n> k = 3mN nt k+i, with N n _ n = 3m. 
Hence N n j. = (3m)™~ fe+1 . It is easily verified that 
Nn = 3 + ELi N n,k and B n = \C n + ± ^Li N n , k D k . 
Note that the number of internal vertices with a given de- 
gree D satisfies N(D) ~ D~ K where k = ln(3m)/ ln(m). 
For the Sierpinski gasket at m = 1, all internal vertices 
have the same coordination number D = 4. 

We now assign a Potts variable <7j = 0, 1, . . . , q — 1 to 
each of the lattice sites. The energy of a Potts configu- 
ration {(Ti} is given by, 



H = -J^% i ,(T J )-JI^%i,0). 

(ij) 1 



(2) 



The first sum extends over all pairs of sites connected by 
a bond, while the second sum introduces a symmetry- 
breaking field that favors or disfavors a selected Potts 
spin state a — 0. Here 8(a, a 1 ) = 1 if cr = a 1 and 
otherwise. In this work we shall restrict ourselves to the 
ferromagnetic case at J > 0. 

Model ([2]) is not invariant under the RG transforma- 
tion discussed below. However, it can be recast in the 
form of an invariant model that contains the one-spin 
term H and a set of three-spin (which includes the pair- 
wise) interactions on a triangle. Making one spin state 
(e.g., o~ = 0) special, there are altogether seven non- 
symmetry-related three-spin configurations as shown in 
Fig. 2. At a given temperature T, the three-spin inter- 
actions can be expressed in terms of Boltzmann weights 
{Z a } for the triplet configurations 12]. The one-spin 
term is represented by a weight v = cxp(H/T) assigned 
to all internal sites in the state a = 0. A restricted 
partition function Z a ^ of the fractal lattice at order fc 
is obtained by summing over spins on all internal sites, 
whereas spins at the three corner sites define the label a 
according to the scheme shown in Fig. 2. 




FIG. 2. (Color online.) Seven non-symmetry-related three- 
spin configurations and the associated Boltzmann weights (or 
restricted partition functions). Open circles represent spins 
in the selected state a = 0. Solid circles represent spins in 
the remaining q — 1 states distinguished by color. 



III. ITERATIVE CALCULATION OF THE 
PARTITION FUNCTION AND THE 
RENORMALIZATION GROUP FLOW 

The partition function of the Potts model with three- 
spin interactions can be computed exactly by following 
the two-step construction of Fig. 1. In the first step, 
{Z' a } on the lattice fc' are obtained from {Z a } by sum- 
ming over spins on the three internal sites (solid circles). 
The mapping {Z a } — > {Z' a } for general q was worked out 
by Borjan et ah, where v enters as a parameter of the RG 
transformation [l2|, EH . In the second step, each of the 
seven weights are replicated as Z' a — > Z a = (Z' a ) m . 

To illustrate the procedure, we consider first the sym- 
metric case where the a = state is made identical to the 
other q - 1 states, i.e., Z xfi = Z 4 , , Z 2 ,o = -^3,0 = ^5,0, 
Z 6 ,o = ^7,0) an d v = 1. This includes model © at H = 
0. The RG transformation defines a two-dimensional 
map for the variables x = Z^jZ± and y = Z7/Z4. Let 
z a = Z' a /Zl, the general iterative relations of Ref. [l2T | 
reduces to: 

z 4 (x, y) = l + (q- l)x 2 (3 + Ax) 

+(q - l)(q - 2)y [3x 2 + 3xy + (q - 3)y 2 ] , 

zs(x, y) = x + Ax 2 + 3x 3 

+(q - 2) (2xy + y 2 + 5x 3 + 8x 2 y + 3xy 2 ) 
+ (q -2){q- 3)y [x 2 + 7xy + (q - <3)y 2 ] , (3) 

z 7 (x, y) = 3x 2 + 6xy + 8x 3 + 9x 2 y + y 3 

+(q - 3) (3y 2 + x 3 + 21x 2 y + 9xy 2 + 3y 3 ) 
+ (q-3)(q-A)y 2 [9x+(q-2)y}. 

The new ratios of the restricted partition functions after 
one iteration are then given by, 



x = 4>{x,y) 
y = i/)(x,y) 



z 5 {x,y)i m 



z 4 (x,y). 

zr{x,y)~ 



(4) 



Specializing on the model at H = 0, we have 
Z 4 o = exp(3mJ/T), Z 5 = exp(mJ/T), and Z7 = 1. 
The starting values of the iteration are xq — ^5,0/^4,0 = 
cxp(— 2mJ/T) and j/o = ^7,0/^4,0 = exp(— 3mJ/T). 
Successive applications of Eqs. <j3j> and (|4|) yield the ra- 
tios Xk+i = 4>{x k ,y k ) and y k +i = ijj(x kl y k ). 

The full partition function on the nth order triangle is 
obtained from the restricted partition functions, 

Z k = q [1 + 3(q - l)x k + (q-l)(q- 2)y k ] Z 4 , k , (5) 

where Z^ k satisfies the iterative relation, 

ZiM = \z±{x k _ u yk-x)Z% k ^\ m . (6) 

To compute the free energy F k = —T\nZ k , we obtain 
from Eq. ©, 

\nZ^k = mhiZ4,(xk-i,yk-i) + 3mlnZ 4ife _i 

fc-i 

— n — 1 

\nz 4 (x n ,y n ). (7) 

The free energy per unit triangle in the infinite size limit 
is thus expressed as, 

p 00 
f(T)= lim -f = -3J-T V(3m)- n - 1 lnz 4 (a; Il ,y„). 

(8) 

The mapping defined by Eqs. ([3]) and ((4]) has two 
trivial fixed points at L = (0,0) and H = (1,1) on the 
(x, y)-plane, corresponding to the ordered and disordered 
phases, respectively. As noted previously (loj. on the Sier- 
pinski gasket at m = 1, the mapping at small x yields, 

xkx + 4:X 2 , (9) 

which grows under the iteration. Therefore the low tem- 
perature fixed point L is unstable and no finite tempera- 
ture transition is expected. 

For m > 1, we have 4>{x,y) ~ x m at small x which 
turns L into a stable fixed point. Consequently, the Potts 
model on such a lattice is ordered at sufficiently low tem- 
peratures 0. A third fixed point C = (x c ,y c ) emerges 
and controls the critical properties at the transition to 
the high temperature disordered phase. Figure 3 shows 
the flow structure of the map at m = 2 and q = 3 which 
is representative of the general situation. The thick solid 
line is the critical manifold that separates the (x, y)-planc 
into two regions. Below this line, the mapping carries the 
system to the low temperature fixed point L. Hence this 
part of the parameter space is associated with the ordered 
phase. Above the line, the mapping carries the system 
to the high temperature fixed point H and the system is 
disordered. On the critical manifold, the mapping carries 
the system to the hyperbolic fixed point C. 

In the neighborhood of the fixed point C, the mapping 
can be well approximated by a linearized version applied 
to 5x k = x k - x c and Sy k = y k - y c , 

( J**-* ) = ( y )(f k ), (10) 

^ 6y k+1 J \ d x tp dyip J \Sy k J y ' 




FIG. 3. Fixed points (solid circles) and invariant manifolds 
(solid lines) of the renormalization group transformation at 
m — 2 and q = 3. The critical manifold (thick solid line) 
divides the parameter space into high and low temperature 
regions. The dashed line given by y = x 3 ^ 2 corresponds to 
the Potts model © at H = 0. 



TABLE I. Fixed points of the RG transformation at v = 1. 
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where the derivatives are evaluated at (x c , y c ). Diagonal- 
ization of the matrix yields two eigenvalues Ai = Xt > 1 
and A2 < 1 that describe the expansion and contraction 
rates along the unstable and stable manifolds of the hy- 
perbolic fixed point C, respectively. 

More generally, using Z4 as the normalization factor, 
the RG transformation applied to the seven restricted 
partition functions yields a nonlinear map in six dimen- 
sions spanned by g a = Z a jZ 4 . At v = 1, the symmetric 
case defines a two-dimensional invariant manifold of the 
mapping, with L, H and C the three fixed points on this 
manifold. By tracing the RG flow along the unstable di- 
rections emanating from these points, we obtain a total 
of 9 fixed points as listed in Table 1. Also given in the 
table is the number of unstable directions n u (known as 
the codimension) and the maximum eigenvalue A m of the 
linearized map at each fixed point. 

Figure 4 illustrates schematically the RG flow between 
the fixed points. Points labelled by L', C, H' and L, 
C", H" lie on the decoupled manifold with 172 = 33 = 
ge = 0. In this case, Z\ = Zf m transforms by itself. 
The remaining three weights Z4, Z5 and Zj transform 
according to Eqs. and ||3J| with the replacement q — > 
q' = q — 1. The critical fixed point of the (q — l)-state 
Potts model is denoted by (x' c ,y' c ). At criticality, Z' 4 
grows by a factor (z' 4 c ) m upon each iteration. 

The above RG flow topology is largely preserved in a 
non-zero field H. The decoupled manifold at gi — 9z — 



#6 = is invariant under the RG transformation at any 
v. Among the two fixed points not on this manifold, H 
is always a fixed point of the transformation. The point 
C survives when H ^ 0, but its location shifts with the 
broken symmetry. 

The structure of the RG flow as shown in Fig. 4 is rem- 
iniscent of the two-dimensional site-diluted three-state 
Potts model for the physisorption of Kr on graphite stud- 
ied by Berker et al. [15[ The three globally stable fixed 
points L', H' and G each represent a "bulk phase": L' 
for the (q — l)-state ordered or in-registry solid phase, H' 
for the (q— l)-state disordered or liquid phase, and G for 
the vacancy (er — 0) dominated gas phase. The shaded 
surface in Fig. 4 indicates the phase boundary between 
the condensed (L' and H') and gaseous (G) phases. This 
transition is first order except on the line CH, to be dis- 
cussed below. On the other hand, transition between the 
ordered and disordered condensed phases (cross-hatched 
surface in Fig. 4), controlled by the (q — l)-state hyper- 
bolic fixed point C', is continuous. 



IV. PHASE DIAGRAM OF THE POTTS 
MODEL AND CRITICAL EXPONENTS 

The phase diagram of the Potts model defined by 
@ follows from the general topological structure of the 
renormalization group flow presented in the previous sec- 
tion. As usual, the nature of ordering in the system at 
a given T and H is determined by which of the three 
fixed points L', H' and G the RG transformation carries 
the set of reduced partition functions to. From symme- 
try considerations, the phase associated with G occupies 
the quadrant at H > 0. On the H < side, an order- 
disorder transition among the q — 1 favored Potts states 
takes place at a temperature T C (H). The transition tem- 



L 




FIG. 4. Topological organization of selected fixed points and 
invariant manifolds under the renormalization group trans- 
formation at v = 1. Solid circles in black are stable fixed 
points associated with the three bulk phases, separated by 
the shaded and cross-hatched surfaces that represent first or- 
der and continuous transitions, respectively. Open circles lie 
on the invariant manifold where all q spin states are identical. 
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FIG. 5. Phase diagram of the Potts model ([2]) at m = 2 and 
q = 3. Lines of first order (solid) and continuous (dashed 
and dotted) transitions meet at the bicritical point H = 0, 
T = T bc . 

perature at a given H can be determined numerically by 
noting that the system flows to L' for T < T C (H) but 
to H' for T > T C (H) under successive RG transforma- 
tions. Note that this transition is absent in the special 
case q = 2. 

Figure 5 shows the phase diagram of the Potts model 
for the representative case q = 3 and m — 2. In the 
absence of the symmetry-breaking field (i.e., H = 0), 
spontaneous ordering occurs for T < T^ c ~ 2.8919 with 
g-fold degeneracy. A positive (negative) H selects (de- 
selects) the a — state. Thus, as in the case of the Ising 
model at q = 2, the line at H = 0, T < T^c (solid line 
in Fig. 5) becomes a first order transition line on the T- 
H plane. The order-disorder transition continues on the 
H < side (dashed line in Fig. 5) and joins the line of 
first order transition tangentially at the bicritical point 
(Tb c ,H = 0)[lj|. Furthermore, the gas and the (q — 1)- 
state disordered phases are separated by a continuous 
transition on the high temperature side. 

We now turn to the calculation of the critical exponents 
that describe the singular behavior of various thermo- 
dynamic quantities at the continuous phase transitions. 
The standard RG procedure allows them to be computed 
from the eigenvalues of the linearized map at the appro- 
priate fixed points. For example, the order-disorder tran- 
sition at H < (dashed line in Fig. 5) is governed by 
the fixed point C which has a single unstable direction 
with an eigenvalue X' T (q) = Xr(q — !)■ In the neighbor- 
hood of the transition point T = T C (H), the free energy 
contains a singular part f s ~ \T — T c (H)\ 2 ~ a . The stan- 
dard scaling analysis yields the specific heat exponent 
a' = 2 — ln(3m)/ In for the (q — l)-state Potts model, 
independent of H. 

The first order transition at H — 0, T < Tb c is gov- 
erned by the fixed point L. The usual scaling argument 
suggests that the singular part of the free energy behaves 
as f s - \H\ 1+1 / 5 , where 1 + 1/6 = lnd//lnA, with A be- 
ing the eigenvalue for the unstable direction of the RG 
flow. At L, A = 3m = df gives 6 = 5^ = oo, hence 

f.~\H\, (T < T bc ). (11) 

Equation can also be derived from the jump in the 
fraction po of vertices in the a = state across the tran- 



sition line. At H = 0, we have "phase-coexistence" of 
the q spontaneous symmetry-breaking states, each favor- 
ing a particular Potts state whose fraction on the lattice 
p(T) > 1/q. A positive field H selects the a — state 
with po — p(T). On the other hand, a negative H de- 
selects the a = state, yielding po = [1— p(T)]/(q— 1). In 
each case, the free energy changes from its H — value 
by an amount A/ = —Hpo per vertex. The jump in po 
gives rise to the singularity in the free energy above. 

For T > T bc , dependence of the free energy on H is 
controlled by the RG flow near the fixed point H. Table 1 
gives A = m and hence 6 = 5n = l/[ln(3m)/ ln(m) — 1] = 
In m/ In 3. Consequently, 

/ s „ | jff |l+(ln3/lnm) ) (T > T hc ). (12) 

Such a singular behavior on the high temperature side is 
quite unusual. The common notion is that the disordered 
system consists of statistically independent subsystems 
of finite size. Since the free energy of each subsystem 
is analytic in the controlling parameters T and H, the 
free energy of the system as a whole should be analytic 
as well. However, examining Eq. ^ more closely, we 
find an interesting possibility. Consider a fcth order ver- 
tex and its CV- = 2m k+1 neighbors. A positive external 
field H > generates an excess probability of order H 
for each neighboring spin to be in the a = state. This 
gives rise to an effective field H e s,k ~ HCk — Hm k at 
the vertex in question from the ferromagnetic couplings, 
which grows exponentially with k. Setting H e g = T 
yields a critical value k c — — hx{H/T)/\nm. Vertices 
with k > k c are completely ordered. Contribution to the 
free energy per triangle from these vertices is estimated 

to be f s ~ Tj^-Heg^ ~ (3m)~ fc - ~ #ln(3m)/lnm which 

agrees with Eq. ([12)) . 

We now consider the most complex situation near the 
bicritical point at H = 0,T = T hc . The RG flow is 
controlled by the fixed point C which has three unstable 
directions as shown in Fig. 4. In terms of the scaling 
fields t = (T — T hc )/T hc along CH, H along CG, and 
He along CC", the singular part of the free energy is 
expected to satisfy the scaling, 

f s (t,H,H c ) = b- d! f s (by T t,by»H,by c H c ), (13) 

where b is an arbitrary scaling factor. The scaling expo- 
nents y a are related to the eigenvalues A a of the linearized 
map through y a = lnA a /ln2. For the model defined by 
Eq. ([2]), He can be expressed as a linear combination of 
t and H but since He = when H = 0, we may write 
H c = H. With the choice b = \t\- x / VT , Eq. (13]) takes a 
more suggestive form, 

m,H) = \tr*r ± (^), d4) 

where the subscript '±' indicates the sign of t. Here 
a = 2 — ln(3m)/ ui(At), A = \nX H /\n\ T , and = 
In Ac/ In At are the specific heat, gap, and crossover ex- 
ponents, respectively. 
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FIG. 6. Eigenvalues \h, At and Ac associated with the three 
unstable directions of the renormalization group flow at C 
against 1/q. Here m = 2. Also shown are the specific heat 
exponent a (dashed line), the gap exponent A (dotted line), 
and the crossover exponent (j> (dash-dotted line). 



Figure 6 shows the numerically computed eigenvalues 
Xh , At, Ac and the exponents a, A, <fi at m = 2 but dif- 
ferent q values. It is evident that A# > At > Ac and 
A > 1 > (f> for all q. Dominance of the symmetry- 
breaking field renders most of the neighborhood of the 
bicritical point to be occupied by either the (q — l)-state 
ordered phase at H < or the gas phase at H > 0. 
The (q — l)-state disordered phase appears in a narrow 
power-law wedge on the H < side, bounded below by 
the transition line H = H C (T) ~ -At A . Within this 
phase, H/t^ < At^-' 1 -' ->• as t -> 0. The scaling func- 
tion J-+(u, w) in Eq. (|14[) is singular on the two transition 
lines at u = and u = —A, respectively. In the latter 
case, the scaled variable w = H/t^ controls the crossover 
from the q-state to the (q — l)-state Potts criticality. De- 
tails will not be elaborated here. 

At large q values, one may approximate the mapping 
near C with 4>(x, y) ~ ip(x, y) ~ (qy) 3m which yields x c ~ 
y c ~ q-3m/(3m-i) _ Consequently, At(<7 ->■ oo) = 3m, in 
agreement with the trend seen in Fig. 6. In this limit, 
C eventually merges with L, yielding Xh — 3m as well. 
Both a and A approach 1. The magnetization exponent 
/3 = 2 — a — A— > 0. We have checked that the above 
behavior is quite general for any m > 1. 



V. SUMMARY 

In this paper, we analyzed in detail solution of the 
Potts model on a fractal lattice that has a low temper- 
ature ordered phase. Compared to the Sierpinski gasket 
at m = 1, the coordination number of vertices at the 
corners of the fractal triangle with m > 1 grows with lat- 
tice size. Thus our findings are consistent with previous 
understanding that a finite T c requires sufficiently strong 
connectivity on the fractal lattice. 

Exact renormalization group transformation of the re- 
stricted partition functions is presented here with the 
help of previous work by Borjan et al. 12] Although 
the mapping in the six-dimensional space contains many 
fixed points and invariant manifolds, a relatively sim- 
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pie RG flow diagram is obtained. The flow is mainly 
denned by three stable fixed points, each representing a 
bulk phase with a distinct symmetry or broken symmetry. 
The topology of the RG flow is robust against the lattice 
geometry parameter m and the number of Potts states 
q. By combining analytical and numerical approaches, 
we have computed the phase diagram of the model and 
a complete set of critical exponents for the transitions, 
including those at the bicritical point, based on the stan- 
dard RG framework. 

It is worth noting that, on the laced Sierpinski gasket, 
transition to the high temperature disordered phase is 
continuous for all q values. Previously, solution of the 
Potts model on the Berker lattice also yielded contin- 



uous transitions with a specific heat exponent a(q) that 
approaches 1 in the limit q —> oo 1 , 7] . This is in contrast 
to the g-state Potts model on periodic lattices in two and 
higher dimensions, where the transition becomes first or- 
der when q exceeds a critical value q c . Further work is 
needed to gain a precise understanding of this difference. 
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